Antiprotonic Potentials from Global Fits to 

the PS209 Data* 



E. Friedman 1,2 , A. Gal 

Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 



Abstract 

The experimental results for strong interaction effects in antiprotonic atoms by 
the PS209 collaboration consist of high quality data for several sequences of iso- 
topes along the periodic table. Global analysis of these data in terms of a p-nucleus 
optical potential achieves good description of the data using a s-wave finite-range 
pN interaction. Equally good fits are also obtained with a poorly-defined zero-range 
potential containing a p-wave term. 
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1 Introduction 



In earlier publications [1,2] most of the then available data on strong interac- 
tion effects in antiprotonic atoms have been used in 'global' analyses in terms 
of a p-nucleus optical potential. The simplest form of the potential is 

2pV opt (r) = -4n(l + ^^-j^)[b (p n + p p ) + b 1 (p n -p p )] , (1) 

where p is the p-nucleus reduced mass, p n and p p are the neutron and proton 
density distributions normalized to the number of neutrons iV and number 
of protons Z, respectively, A = N + Z and M is the mass of the nucleon. 
In the zero-range Hp* approach the parameters &o and b\ are minus the p- 
nucleon isoscalar and isovector scattering lengths, respectively, otherwise these 
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parameters may be regarded as 'effective'. The fits of Ref. [1] were based on 
48 data points which included only two pairs of isotopes. As a typical result 
we quote Reb = 2.51 ± 0.24 fm, lmb = 3.46 ± 0.27 fm, leading to x 2 / N > the 
X 2 per point, of 1.5, which represents a very good fit to the data. Additional 
terms had also been cosidered in [1] such as a p-wave term, which led to x 2 /N 
in the range of 1.2 to 1.4. Such a reduction in x 2 is, however, not statistically 
significant in these circumstances. Consequently only the above zero-range bo 
version has since been used in several analyses of experimental results. 

The recent publication of results for strong interaction effects in antiprotonic 
atoms by the PS209 collaboration [3] changed the experimental situation dra- 
matically by providing high quality data for several sequences of isotopes along 
the periodic table. Such data hold promise to pin-point isospin effects in the 
potential and/or explicit surface effects such as the existence of a p-wave term 
in the potential. This is the scope of the present work which is based on data 
for five isotopes of Ca, four isotopes each of Fe and of Ni, two isotopes each of 
Zr and of Cd, four isotopes of Sn, five isotopes of Te and the isotope 208 Pb, all 
from results of the PS209 collaboration [3]. Note that the results for Cd, Sn 
and Te isotopes have been revised very recently to include effects of E2 reso- 
nance [4,5]. In order to expand the range of the data further we have included 
also earlier results for 16 > 18 (see [1] for references to the original experiments.) 
In section 2 we discuss the choice of nuclear density distributions p n and p p , 
which form an essential ingredient of the optical potential Eq. (1). Section 3 
summarizes the results and presents a 'recommended' set of parameters which 
best describes the data. 



2 Density distributions 

In choosing the nuclear density distributions for use in the analysis of strong 
interaction effects in antiprotonic atoms, one is faced with the fact that due 
to the strength of the p-nucleus interaction and the dominant role played by 
the imaginary potential, strong interaction effects reflect the interaction with 
the extreme outer regions of the nucleus where the density is well below 10% 
of the central density. It is therefore not very reliable to infer from these data 
quantities which relate to the bulk of the density such as root mean square 
(rms) radii. In what follows we therefore rely also on information from pionic 
atoms [6] when choosing the densities, bacause pionic atoms are sensitive to 
the region where the density is 50% or more of the central density. The density 
distribution of the protons is usually considered known as it is obtained from 
the nuclear charge distribution [7] by unfolding the finite size of the charge of 
the proton. The density distribution of the neutrons is usually not known to 
sufficient accuracy. 
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Fig. 1. Four neutron and two proton densities for 120 Sn. Solid curves for SP model, 
dotted curves for 'skin' form, dashed curve for 'halo' form and dot-dashed curve for 
the 'average' form. All neutron densities have the same rms radius. 

In a recent paper [8] Trzciriska et al. have suggested, based on analyses of 
strong interaction effects in antiprotonic atoms and of radiochemical data 
following p absorption on nuclei, that the difference between neutron and 
proton rms radii in nuclei depend on the asymmetry parameter (N — Z)/A as 
follows: 

N — Z 

r n -r p = -(0.04 ±0.03) + (1.01 ±0.15) — - — (fm). (2) 



Addressing the question of the shape of the neutron density, they have con- 
cluded that a 'halo' rather than a 'skin' shape is preferred. These shapes are 
defined as follows: Assuming a two-parameter Fermi distribution both for the 
proton (unfolded from the charge distribution) and for the neutron density 
distributions 

PnAr) = l + eM(r-Rn, P )/a n , p y (3) 



then for the 'skin' form the same diffuseness parameter for the protons and the 
neutrons a n = a p is used and the R n parameter is determined from the rms 
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radius r n . In the 'halo' form the same radius parameter R n = R p is assumed for 
p n and for p p , and is determined from r n . We also use an 'average' option 
where the diffuseness parameter is set to be the average of the above two 
diffuseness parameters, namely, a^ ve = (a p + a^)/2 and the radius parameter 
R n is determined from the rms radius r n . To allow for different values of r n — r p 
we have used Eq. (2) but with variable slope parameter, where, e.g. a slope of 
1.5 fm correspons to the average behaviour of results of relativistic mean field 
(RMF) calculations [6,9]. 

Figure 1 shows examples of densities used in the present work. The single- 
particle (SP) densities are described in [2]. The various neutron densities have 
the same value of r n , as given by RMF calculations. Note that where the 
densities are well below 10% of their central values the neutron densities are 
significantly larger than the proton density. 



3 Results 

Fits were made to 90 data points using a variety of shapes of nuclear densities 
and values of r n as described in the previous section. Finite range was also 
introduced as an option, where the nuclear densities have been folded with a 
Gaussian, separately for the real and for the imaginary potentials, such that 
the densities in Eq. (1) are replaced by 'folded' densities 

f(r) = / rfrVM^e-C-'W . (4) 

Table 1 summarizes the main features of the results. The columns marked 
'RMF' were obtained with values of r n — r p taken from results of RMF cal- 
culation [9] whereas the columns marked £ p' are for r n — r p due to Ref. [8]. 
In both cases the 2-parameter Fermi distribution was used for the proton and 
for the neutron densities. The column marked with SP is for both densities 
obtained by filling in levels in single particle potentials [2] with r n — r p values 
from the RMF model. The shapes of the neutron distributions ('ave.', 'halo', 
'skin') are as explained in the previous section. The top part of the table is 
for zero range interaction and the bottom part of the table is for finite range 
interaction with Gaussian range parameters /3r and (3j for the real and for 
the imaginary parts, respectively. From the values of x 2 /N for the zero range 
fits it is seen that reasonable fits to the data can be achieved for all options 
except for the 'RMF-halo' variety. The increased values of x 2 /N compared to 
results of analyses of the previous data [1,2] reflect the improved accuracy and 
the additional information contained in the new data. The latter is observed 
in the lower part of the table where the need for a finite-range version of the 
p-nucleus optical potential is self-evident. The values of x 2 and of x 2 /N mean 
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Table 1 

Global fits to 90 antiprotonic atom points. 'RMF' and 'p' stand for values of r n — r p 
used. Top part is for zero-range potentials and bottom part for finite-range potentials 
with range parameters (3rj- 

'RMF' 'p atoms' 





SP 


ave. 


halo 


skin 


ave. 


halo 


skin 


x 2 


296 


338 


689 


242 


240 


310 


249 


X 2 /N 


3.3 


3.8 


7.7 


2.7 


2.7 


3.5 


2.8 


Re6o (fm) 


2.5(2) 


2.6(2) 


2.0(2) 


3.4(2) 


3.1(2) 


2.7(2) 


3.6(2) 


Im6o (fm) 


3.1(2) 


2.8(2) 


2.5(2) 


3.1(2) 


3.2(2) 


3.1(2) 


3.4(2) 


x 2 


221 


203 


209 


219 


206 


202 


235 


X 2 /N 


2.5 


2.3 


2.3 


2.4 


2.3 


2.2 


2.6 


Re6o (fm) 


0.78(5) 


1.20(7) 


0.68(5) 


1.11(5) 


1.77(9) 


1.34(7) 


0.71(3) 


Im^o (fm) 


1.75(6) 


1.74(8) 


1.10(4) 


2.8(1) 


2.3(1) 


2.0(1) 


3.5(1) 


Pr (fm) 


1.0 


0.9 


0.9 


1.2 


0.7 


0.8 


1.6 


Pi (fm) 


1.2 


0.9 


1.3 


0.6 


0.7 


0.8 


0.6 



significant improvements in the fits compared to the zero-range results. The 
best values for x 2 /N are close to what has been achieved in global fits to pionic 
atom data [6] . 

Next we turn to the question of whether the data require the presence of an 
isovector term b\ in the potential Eq. (1). Fits to the full data set of 90 points 
failed to produce any reduction in x 2 when such a term was introduced. In 
a separate calculation we have been able to obtain a non-zero value for the 
parameter b\ when fits were made to only the data for the Sn isotopes, leading 
to % 2 /N of 1.0 for these isotopes. However, applying that b\ to the full data set 
resulted in a major deterioration of the fit, and a further global fit restored the 
results of Table 1 with &i=0. This is not surprizing since antiprotonic atoms 
are sensitive to the nuclear density in the extreme outer regions of the nucleus 
where the neutron density is dominant. Explicit isospin effects are, therefore, 
hard to observe. 

Another attempt to improve the agreement between zero-range calculations 
and experiment was made by introducing a p-wave term into the optical po- 
tential [1,2], a term which is active mostly at the surface region. Starting with 
the zero-range version of the potential, the values of \ 2 could be reduced from 
those of the upper part of Table 1 to those of the lower part, but with p-wave 
parameters which depended critically on the choice of shape for the neutron 



5 



distribution. Including in addition finite-range folding failed to further improve 
the fits. 



In conclusion, global analyses of most of the data from PS209 reveal a need 
to improve on the previously successful zero-range s-wave isoscalar p-nucleus 
optical potential. Including an isovector term is not required by the data nor 
is a p-wave term recommended at our present state of knowledge, although 
it undoubtedly improves the fits to the data when the zero-range form is 
used. The use of a finite-range interaction in the optical potential is definitely 
required by the data. Assuming that the values of the differences between 
rms radii r n — r p are, on the average, between those advocated in Ref. [8] 
and those predicted by RMF calculations [9], and assuming the shape of the 
neutron distributions is between those described by the 'skin' and by the 
'average' prescriptions, we then recommend the values Re&o = 1.3±0.1 fm, 
Im&o = 1.9±0.1 fm, with /3r = Pi = 0.85 fm within the folding expression Eq. 
(4). 
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